Magnetoelectric Response of Multiferroic BiFeC>3 and Related Materials 
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We present a first-principles scheme for computing the magnetoelectric response of multiferroics. 
We apply our method to BiFe03 (BFO) and related compounds in which Fe is substituted by other 
magnetic species. We show that under certain relevant conditions - i.e., in absence of incommensu- 
rate spin modulation, as in BFO thin films and some BFO-based solid solutions - these materials 
display a large linear magnetoelectric response. Our calculations reveal the atomistic origin of the 
coupling and allow us to identify the most promising strategies to enhance it. 

PACS numbers: 75.80.+q, 71.15.Mb 

Magnetoelectric (ME) multiferroics present coupled 
electric and magnetic orders jlj, which may allow the 
development of a variety of magnetic devices, from mem- 
ories to spin filters, whose behavior would be switchable 
by application of a voltage. The experimental search for 
robust room-temperature (T room ) multiferroics is prov- 
ing a major challenge. The first-principles contribution 
to this effort is quickly increasing, as shown by recent pre- 
dictions of new materials @, 0, |j| and novel ME coupling 
mechanisms 0, @, 0]- Yet, many key issues remain to be 
addressed theoretically, as even the ME response of the 
most promising systems (e.g., thin films of BiFe03 Q) 
still needs to be characterized and understood in detail. 
Here we describe a method for ab initio computations 
of the ME response of multiferroics, tackling the strain- 
mediated contributions that typically exist in these ma- 
terials. We apply our method to BiFeC>3 and related 
compounds. 

Methodology - The response of a linear magnetoelec- 
tric that is not multiferroic - and is thus paraelectric - 
can be split in two contributions a purely electronic 
one and an ionic one. The ionic part accounts for the 
structural response to an electric (magnetic) field and the 
resulting change of magnetization (polarization). When 
the linear magnetoelectric is multiferroic - and thus fer- 
roelectric - there is a third contribution associated to 
the piezoelectric response to the electric field. Indeed, 
it can be shown that all commensurate ll| multiferroics 
must be piezomagnetic, as the symmetries that preclude 
piezomagnetism are broken in such systems [lij ]. Hence, 
the combination of piezoelectricity and piezomagnetism 
in multiferroics results in a strain-mediated contribution 
to the linear ME response. 

In order to model such effects, we have generalized to 
the magnetoelectric case the formalism that Wu, Van- 
derbilt, and Hamann [l3| (WVH in the following) intro- 
duced for a systematic ab initio treatment of dielectric 
and piezoelectric responses. Let us consider how the en- 
ergy of a multiferroic crystal varies as a function of the 
following perturbations: displacements u m of atoms away 
from their equilibrium positions, homogeneous strains rjj, 
and applied electric (£ a ) and magnetic (TL^) fields. Here, 
to is a composite label that runs over the atoms in the 
unit cell and displacement directions, j runs from 1 to 6 
in Voigt notation, and a and fi are Cartesian directions. 



We write the energy per undeformed unit cell volume of 
the multiferroic crystal, expanded to second order around 
its zero-field equilibrium structure, as 
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assuming summation over repeated indexes. The equi- 
librium condition implies the terms linear in u m and rjj 
are zero; the term A a ——P^ [A^—^M.^) is the sponta- 
neous polarization (magnetization). The second deriva- 
tives of Eq. ([1]) correspond to well-known physical quan- 
tities, such as the force constant matrix K mn — QoB mn , 
the frozen- ion elastic tensor Cjk = Bjk, the Born dynam- 
ical effective charges Z ma = —floB moci the force-response 
internal strain tensor A mj = —£loB m j, or the frozen-ion 
piezoelectric stress tensor e QJ = —B a j, where f2 is the 
unit cell volume and the bar stands for frozen-ion. Ad- 
ditionally, the purely electronic ME response is 
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where the hat stands for frozen-cell. 
piezomagnetic stress tensor is 



The frozen-ion 
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Lastly, the magnetization change driven by an atomic 
displacement is 
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The tensors describing the response to static fields should 
account for the field-induced ionic relaxation. To com- 
pute such relaxed-ion quantities, we introduce the energy 
E(£,Tt,ri) = min u E(u, £, Tt, rf}, which is obtained from 
Eq. Q]) by replacing 



(B) mn (B na £ a 



B 



njVjJ 



(5) 



E can be written in a form analogous to Eq. (p} , with B 
coefficients that are combinations of the original £?'s. By 
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computing the second derivatives of E we can obtain, for 
example, the piezomagnetic stress tensor 



h u 



d 2 E 



-B u 



or the frozen-cell ME tensor 
d 2 E 



(6) 
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Finally, in order to compute the full ME response we 
go one step beyond WVH and introduce the energy 
£(£ , 7i) = min,, , Ti, rf), from which we derive 
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Note that, except for the purely electronic a™, all the 
terms in Eq. © can be obtained from finite difference 
calculations that do not require the simulation of applied 
fields. Thus, the present formalism brings the study of 
ME effects in multiferroics within the scope of the most 
widely used Density Functional Theory (DFT) codes. 

Finally, let us remark we would like our theory to de- 
scribe the so-called proper piezoelectric effects. WVH 
dealt with this issue by introducing suitably rescaled elec- 
tric fields and polarizations [Hj]. Equivalently, the de- 
sired result is obtained if the B a j coefficients of Eq. (UJ 
are computed by finite differences as described in Ref. 1 14 . 
Interestingly, the proper vs improper distinction applies 
to the piezomagnetic case too, a fact that has seemingly 
remained unnoticed. Indeed, we want our proper B^j co- 
efficients to be such that (i) H-induced rotations of the 
sample do not contribute to the response and (ii) a mere 
change in the unit cell volume has no magnetic effect. 
The former requirement is easily accomplished by work- 
ing only with the symmetric part of the strain tensor, 
and we comply with the latter by computing the B^j 
coefficients from strain-induced changes in the magnetic 
moment per cell, as opposed to changes in the magneti- 
zation [151 ]. 

Application to BiFeO^ and related materials - We ap- 
plied our method to BiFeC>3 (BFO), arguably the most 
promising multiferroic. We simulated BFO with the mag- 
netic structure of the systems that are most relevant for 
applications, i.e., BFO thin films and solid solutions in 
which Bi is partly substituted by a lanthanide to im- 
prove the dielectric properties. In such cases BFO looses 
the spin cycloid that occurs in bulk samples (lil . H?! Il8j 
and presents a canted G-type antiferromagnetic (AFM) 
spin arrangement. The magnetic easy axis lies within 
the plane perpendicular to the polar direction (xy in 
Fig. [IJa), z being parallel to V s and corresponding to 
the [llljpc pseudocubic direction) [It], [HJ, and the cor- 
responding magnetic space group allows for a linear ME 
response [20|. Indeed, there is experimental evidence of 



TABLE I: Top: Frozen-cell and full ME tensors for three 
representative cases. Results given in 10 -4 Gaussian units 
(g.u.). Magnetic space groups and easy axes are indicated. 
Bi2FeCr06 has nearly R3 symmetry, except for a tiny spin 
canting. Dots indicate coefficients that are zero by symmetry. 
Bottom: Maximum ME response, as quantified by square root 
of largest eigenvalue of a'ct, for all considered substitutions. 
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such a linear ME effect 21, 22], but we still lack a de- 



tailed and well-established characterization. 

We also studied situations in which, maintaining the 
10-atom cell and basic structure of BFO, Fe is substituted 
by other magnetic species. In particular, we have consid- 
ered Fe^Cr, Fe— >Mn, and Fe^Co substitutions, as well 
as Fc-^(Fe,Cr) and Fe— >(Mn,Ni) double perovskites [24 1. 
Note that, except for Bi 2 FeCr0 6 , none of the studied sub- 
stitutions corresponds to thermodymically stable phases. 
The purpose of considering such ficticious compounds 
(denoted with italics) was to identify chemical trends in 
the magnitude of the ME effects. 

For the calculations we used the Projector Augmented 
Wave approach to DFT, within the so-called "LDA+C/" 
approximation, as implemented in the Vienna Ab initio 
Simulation Package [25]. For all the compositions stud- 
ied, we proceeded as follows: Taking the usual atomic 
and magnetic structure of rhombohedral BFO as start- 
ing point, we relaxed the system to find a well-defined 
energy minimum. This energy minimization included a 
careful search for the magnetic easy axis. Having identi- 
fied the equilibrium state, we computed the coefficients 
in Eq. ([TJ) by finite differences. All the calculations were 
fully self-consistent, allowed for non-collinear magnetism, 
and included spin-orbit couplings. We neglected orbital 
magnetization as well as the purely electronic ME re- 
sponse a, which is likely to be relatively small [2^|, spe- 
cially for materials with a significant structural response 
to electric fields. 

Table Q] summarizes our results and Fig. [lja) sketches 
the f-induced spin canting that underlies the computed 
effects. The obtained ME tensors satisfy symmetry 
relationships within the calculation accuracy; e.g., for 
BiCrOs we got a\2 = —0121, as expected for the 3m 
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magnetic point group [10j . Interestingly, the largest 
ME coefficients obtained share the following two fea- 
tures: (i) The applied £ is perpendicular to V s , which re- 
flects the easy-polarization-rotation mechanism that de- 
termines the largest electromechanical responses in fer- 
roelectrics [27|. (ii) The change in M is perpendicular 



to the direction of the easy axis, which is a signature of 
spin canting and an underlying spin-orbit mechanism. In 
contrast, many symmetry- allowed coefficients were com- 
puted to be essentially null, which reflects (i) the struc- 
tural hardness of specific directions and (ii) the energy 
cost of changing the magnetic moment along the direc- 
tion of the easy axis (which would typically require charge 
transfer between magnetic ions Q). 

For all compounds, the full ME response is bigger than 
the frozen-cell effect, the increase being significant in 
BiCr0 3 , BiMn0 3 , and BiCo0 3 . To illustrate how the 
strain-mediated response operates, we consider the rep- 
resentative case of the ctu coefficient of BiCr0 3 . There, 
the largest strain response to £ i is associated with the 
shear 775, as quantified by the piezoelectric strain coeffi- 
cient G?i5 = eij(C )js — 15 pC/N. The resulting change 
in M2 is quantified by the piezomagnetic stress coeffi- 
cient h,25=— 5xl0 4 A/m. This leads to a strain-mediated 
contribution to the ME response of about — 3xl0~ 4 g.u. 

Our calculations allow us to identify the atomic relax- 
ation mechanisms that mediate the ME response. Let us 
consider first the frozen-cell ME tensor (Eq. ([?]))■ whose 
lattice part can be expressed in terms of the eigenmodes 
of K as ^0 1 Psa K J 1 P^ 1 - Here, k s is the mode eigenvalue 
or stiffness, and p^ a and pf are, respectively, the dielec- 
tric and magnetic polarities introduced in Ref. 0. This 
expression allows us to make a cumulative plot as the 
one shown in Fig. QJb) for S23 of BFO. Clearly, only 
two sets of modes contribute significantly: (i) low-lying 
Bi-dominated modes with a relatively small magnetic po- 
larity (~78 % of the response) and (ii) high-energy Fe- 
dominated modes that involve a relatively large change 
in M. We can thus conclude that the ME response is 
largely driven by Bi modes that do not involve the mag- 
netic species significantly, a result that applies to all the 
studied compounds. The role of Bi is also preeminent in 
the strain-mediated part of the response, the largest B m j 
couplings being unequivocally associated with Bi. 

The work of Ref. [28| allows us to compare theory and 
experiment for BFO. These authors studied the bulk 
material under magnetic fields up to 25 T, which al- 
lowed them to eliminate the spin cycloid and reach a 
phase (G-type AFM) that resembles the one we have 
simulated. One can thus extrapolate the M.(Ti) curve 
measured within the high-field phase back to 7i=0, and 
estimate the magnetization at zero field. The result is 
0.25 emu/g = 0.028 /ie/cell, in remarkable agreement 
with our computed 0.036 /is/cell. Further, the ME re- 
sponse of the high-field phase is about 4xl0 -4 g.u. at 
10 K, again in excellent agreement with our computed 
5xl0 -4 g.u. While this comparison is not fully justified, 
it certainly supports the physical soundness of our rc- 




FIG. 1: Panel (a): Sketch of the response associated to 023 
of BFO. Only the two Fe atoms in the usual (rhombohedral) 
cell of BFO are shown. In equilibrium (dashed arrows) the 
spin of the Fe ions lies on the xy plane; the easy axis is x 
and a small canting results in a net My. When £ y is applied 
we obtain an additional canting and a non-zero M z (solid 
arrows). Panel (b): Cumulative plot for the Q23 coefficient of 
BFO (see text). Modes dominated by Bi (green diamonds) 
and Fe (red squares) displacements are highlighted. 



suits. Comparing our theory with measurements of BFO 
thin films |21| is unfortunately not possible, as the stud- 
ied systems display effects (e.g., very large net magnetic 
moments and possible presence of Fe 2+ ) that clearly do 
not correspond to our simulations. 

To put our results in perspective, note that the largest 
ME responses measured for transition-metal compounds 
correspond to boracites, with a max =20x 10 -4 g.u. for 
Co 3 B 7 Oi 3 Br (see Table Table 1.5.8.2 of Ref. [rJ. Larger 
responses are observed in rare earth compounds, the 
greatest one being a max =100x 10~ 4 g.u. for TbP04. It 
is important to realize, though, that these maximum ct's 
correspond to (very low) temperatures slightly below the 
magnetic ordering transition (17 K for Co3ByOi3Br and 
2.2 K for TbP0 4 0), where the ME effect is strongly 
enhanced. In contrast, the response of BFO computed 
at K is of the same order of magnitude as the a max 's 
of boracites, is expected to grow with T [2^], and should 
occur at T room and above. Hence, when compared with 
other magnetoelectrics, BFO is clearly an unique, very 
promising material. 

Interestingly, from a materials-design perspective one 
could say that BFO is rather unsatisfactory, for several 
reasons. (1) The distortions mediating the ME effect are 
dominated by non-magnetic ions, which is obviously not 
ideal. (2) These materials are quite stiff, as reflected 
by the relatively small piezoelectric and lattice-mediated 
dielectric responses: The largest ones obtained were, re- 
spectively, d max «60 pC/N for BiMn0 3 and elf a t x «45 for 
BiCo0 3 , while for prototype ferroelectric BaTiOa one 
gets d max «250 pC/N and e^R^O [H. (3) These ma- 
terials are poor piezomagnets: For the so-called piezo- 
magnetic strain tensor g l _ l j=h IJ ,k(C~ 1 )jk we obtained a 
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maximum value of 0.4xl0 -10 Oe^ 1 for BiCrOs, which is 
about one order of magnitude smaller than what is typ- 
ical in transition-metal compounds (see Table 1.5.7.2 of 
Ref. [Toh . Yet, in spite of all these drawbacks, BFO is a 
competitive magnetoelectric. 

If we were able to improve on some of these aspects, 
BFO might become an excellent magnetoelectric. Our re- 
sults indicate the ME response of the studied compounds 
owes its relatively large value to the structural response 
to electric fields. The secondary role of the magnetic ef- 
fects is clearly reflected in the fact that the magnitude 
of a max 's (Table QJ does not correlate with the magni- 
tude of the spin-orbit coupling of the considered mag- 
netic species, which rules out chemical substitution of 
iron as a direct way to enhance the ME response. It 
seems more promising to try to increase BFO's electrome- 
chanical responses, which might be achieved, for exam- 
ple, in (i) BFO thin films strain-engineered to be mon- 
oclinic (as opposed to rhombohedral), since monoclinic- 
ity is usually accompanied by structural softness [27j , or 
(ii) Bii-^La^FeOs solid solutions where structural tran- 
sitions occur for small La concentrations [lij]- We hope 
our work will stimulate detailed experimental studies of 
these and similar systems, which our theory suggests con- 
stitute the best hope for obtaining very large linear ME 
effects in BFO. 
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